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Abstract 



This thesis is contributed to the study of decoherence dynamics of the 
dissipative qubit system. We mainly concentrate on the profound impact of 
the formation of a bound state between the qubit and its local environment 
on the decoherence behavior of the reduced qubit system under the non- 
Markovian dynamics. 

Firstly, we evaluate exactly the non-Markovian effect on the decoherence 
dynamics of a single qubit interacting with a dissipative vacuum reservoir. 
We find that the quantum coherence of the qubit can be partially trapped in 
the steady state when the non-Markovian memory effect of the reservoir is 
taken into account. Our analysis shows that it is the formation of a bound 
state between the qubit and its reservoir that results in this residual coherence 
in the steady state under the non-Markovian dynamics. A physical condition 
for the formation of the bound state is given explicitly. Our results suggest 
a potential way to decoherence control by modifying the system-reservoir 
interaction and the spectrum of the reservoir to the non-Markovian regime 
in the scenario of reservoir engineering. 

Secondly, We study the entanglement dynamics of two qubits locally in- 
teracting with their reservoirs and explore the entanglement preservation 
under the non-Markovian dynamics. We show that the existence of a bound 
state of the qubit and its reservoir and the non-Markovian effect are two 
essential ingredients and their interplay plays a crucial role to preserve the 
entanglement in the steady state. When the non-Markovian effect is ne- 
glected, the entanglement sudden death is reproduced. On the other hand, 
when the non-Markovian is significantly strong but the bound state is ab- 
sent, the phenomenon of the entanglement sudden death and its revival is 
recovered. Our formulation presents for the first time a unified picture about 
the entanglement preservation and provides a clear clue on how to preserve 
the entanglement in quantum information processing. 

Finally, in order to obtain a thorough understanding of the entanglement 
dynamics, we study the entanglement distribution of a two-qubit system, 
each of which is embedded into its local reservoir, among all the bipartite 
subsystems including qubit-qubit, qub it-reservoir, and reservoir-reservoir. 
Different to the result that the entanglement of the qubits is transferred 
entirely to the reservoirs under the Markovian dynamics, we find that the 
entanglement can be stably distributed among all components under the 



noil- Mar kovian dynamics, and particularly it also satisfies an identity firstly 
given by Yonac, Yu and Eberly [J. Phys. B 40, S45 (2007)] for a double 
J-C model without decoherence. While the explicit distribution of the en- 
tanglement is dependent of the detail of the model, even the approximation 
used, the identity remains unchanged. Our unified treatment includes the 
previous results in the literature as special cases. The result reveals the pro- 
found nature of the entanglement and should have significant implications 
for quantum information processing. 

This thesis may give a clear clue of decoherence dynamics under different 
approximations and how to preserve quantum coherence in the steady state. 



Contents 



Decoherence dynamics of a dissipative qubit 

2.1 Introduction 

2.2 Tlie model and exact decolierence dynamics of tlie qubit . . 

2.3 Purity and decoherence factor 

2.4 Numerical results and analysis 

2.4.1 The influence of coupling constant 

2.4.2 The influence of cutoff frequency 

2.4.3 The physical mechanism of the decoherence inhibition: 
the formation of atom-photon bound state 

2.5 Summary 



Abstract 

1 Introduction [l| 



12 
U 
U 
16 



11 
20 



Mechanism of entanglement preservation |22 

3.1 Introduction 22 

3.2 The model and entanglement dynamics 24 

3.3 Mechanism of entanglement preservation 26 

3.4 Summary |33 



Entanglement distribution and its invariance 1^5 

4.1 Introduction 

4.2 The model of two qubits in two uncorrelated band-gap reservoirs 

4.3 Entanglement distribution among bipartite subsystems .... 

4.4 Summary 



35 
36 

38 

44 



Summary and outlook |45 



Bibliography |47 

Publication list |52 

Acknowledgments |53 



Chapter 1 
Introduction 



The superposition rule of quantum state, one of the fundamental principles of 
quantum mechanics, allows a quantum system to be in a linear and coherent 
superposition of all possible quantum states pj. It leads to the quantum 
coherence, which is the essential difference of quantum world to the classical 
one. Quantum coherence is of great importance not only in understanding the 
basic rules of quantum mechanics but also in quantum information science 

i- 

Entanglement (also named as quantum correlation), as a non-local quan- 
tum coherence, is one of the characteristic trait of quantum mechanics [T]. 
A state in multipartite system is entangled when it cannot be written as the 
summation of the product states of the subsystems. On the one hand, en- 
tanglement relates to a lot of fundamental problems in quantum mechanics, 
such as, reality, local realism, hidden variables, and quantum measurement 
theory [31 HJ [5]. On the other hand, entanglement can be used as a kind 
of information resource to realize some missions of quantum information 
processing which are intractable for classical one, such as quantum commu- 
nication [6l[7l|8], quantum computation [9l[10], and quantum cryptography 

Any realistic quantum system inevitably interacts with its surrounding 
environment which leads to the loss of quantum coherence, or decoherence of 
the quantum system. This ubiquitous phenomenon deteriorates the super- 
position and the entanglement of quantum state. In terms of information or 
energy, decoherence means that the information or energy flows from quan- 



turn system to the environment irreversibly. The decoherence is deemed as 
one of the main obstacles to the realization of quantum information process- 
ing. Recently, much attention has been paid on the study of the dynamics of 
open quantum system, by which people want to get a thorough understand- 
ing to the detrimental effect caused by decoherence on quantum information 
processing and some clues on how to suppress this unwanted effect. In this 
thesis, we will concentrate on a detailed study of the decoherence dynamics of 
qubits influenced by their vacuum reservoirs under different environments or 
approximations and explore the potential dynamical suppression mechanism 
to the decoherence. 

The dynamics of open quantum system is rather complicated because of 
the complex structure of the environment with which the system of interest 
interacts. Actually, the exactly solvable models are very few, only includ- 
ing the quantum Brownian motion and the system of a two-level atom in 
a vacuum reservoir with Lorentzian spectrum density [Ij. Many approxi- 
mations are ordinarily performed. A generally used approximation in the 
conventional investigation to the dynamics of open quantum system is the 
Born-Markovian approximation [12], which treats the interaction between the 
quantum system of interest and its environment perturbatively and neglects 
the memory effect of the environment. This approximation is valid when 
the coupling between quantum system and its environment is weak (Born 
approximation) and the environmental correlation time is small compared to 
the typical time scale of the quantum system (Markovian approximation). 
This approximation yields equations of motion such as Redfield or master 
equation, which is local in time and mathematically tractable, for the quan- 
tum system of interest. Based on this approximation, it is widely accepted 
that the quantum coherence of a single quantum system flows irreversibly 
to the environment and the decoherence dynamics can be simply depicted 
as an exponential decay. However, things are changed dramatically when 
entanglement dynamics which involves more subsystems, such as qubits, in 
the quantum system, are studied, some works have showed that the entan- 
glement between two qubits ceases abruptly in a finite time scale [131 E]- 
This remarkable phenomenon that the entanglement of the qubits under the 



Markovian decoherence dynamics can be terminated in a finite time despite 
the coherence of single qubit lossing in an asymptotical manner is named 
as entanglement sudden death (ESD). Further investigation shows that such 
ESD is strongly related to the initial portion of double excitation component 
[15J. The larger the initial portion of the component is, the shorter the death 
time is. On the other hand if the environments are composed of the thermal 
or squeezed reservoirs, it is found that the ESD would always happen for any 
initial entangled state. Experimentally, the ESD has been observed using an 
all optical setup and atomic ensemble system [161 IH] • 

The Born- Markovian approximation simplifies greatly the mathematics 
to solve the dynamics of open quantum system, but it suffers more and more 
challenges under the newly emerging experimental results [HI [191 1201 [211 [22] . 
Especially, when the environment has certain structures, such as atom in 
cavity or photonic band gap (PBG) mediums [23l [Ml ESI [261 [271 [28], the 
non-Markovian effect can not be neglected anymore. The non-Markovian 
effect is a kind of dynamic feedback effect which arises from the memory 
effect of environment. In terms of information or energy, the non-Markovian 
effect means that the information or energy flows back from the environment 
to the quantum system of interest. The study of open quantum dynamics 
beyond the Markovian approximation is rather complicated, which needs the 
solving of coupled integro-differential equations. The decoherence dynam- 
ics of quantum system in this case exhibits a dramatic deviation from the 
exponential decay behavior. Actually any kind of environment should have 
memory effect. When this memory effect is very weak, the Markovian ap- 
proximation is applicable. On the other hand, when the memory effect of 
the environment is extremely strong, it would partially feed the lost coher- 
ence back to the quantum system. In this case the Markovian approximation 
would be not applicable. As far as the entanglement in two-qubit system is 
concerned, the non-Markovian effect also has a great impact on it. Model- 
ing the environments as vacuum lossy cavities, Bellomo et al. showed that 
the entanglement would revive again after a finite period time of completely 
disappearance [29]. This is a solvable model and the entanglement dynam- 
ics can be analytically expressed. Via tuning relevant parameters, one can 



easily observe that the non-Markovian effect postpones greatly the death of 
the qubits entanglement. Entanglement dynamics of continues variable sys- 
tem has also been well studied and the non-Markovian effect also makes the 
entanglement dynamics oscillate [301 EB ESj [33l |3l], which can be under- 
stood as the backaction effect of the environment on the quantum system. 
The ESD and its revival due to the non-Markovian effect has been experi- 
mentally observed [22]. All these experimental and theoretically works show 
clearly that the coherence or entanglement time of the quantum system can 
be much enhanced by the non-Markovian effect. 

However, in many cases such finite extension of the coherence/entanglement 
time is not enough for quantum information processing and thus it is desired 
to preserve a significant of the quantum coherence/entanglement, even par- 
tially, in the long time limit forever. Actually, some work has shown that 
it is realizable for some special environment cases. It has been found that 
the spontaneous emission of a two-level atom can be inhibited and its quan- 
tum coherence can be preserved when the atom is placed in a PEG material 
[26| [35| [36] . In the PEG material, the photonic mode density is zero within 
the PEG and this would be accompanied by the classical light localization 
and a photon-atom bound state. The excited-state population in this case 
is partially trapped, a phenomenon known as population trapping [27]. This 
result has been verified experimentally for quantum dot embedded in PEG 
material [37]. It has been reahzed that trapping the single-qubit population 
is the key step to protect entanglement in two-atom case [28]. When two 
initially entangled qubits are immersed in two separate PEG mediums, the 
entanglement of the two qubits can be preserved in the steady state with a 
large fraction. However the mechanism of entanglement preservation is still 
unclear. Also is this a general phenomenon in open quantum system or only 
available in this specific structured environment still is an open question. 

To explore these questions, one should know the dynamics of quantum 
system not only in the short-time scale, but also in the long-time situation. In 
the short-time scale, when the non-Markovian effect is very strong, the quan- 
tum coherence would surfer transient oscillations manifesting the backaction 
effect of the memory environment. It is just the counteraction role played 



by this backaction effect to the dissipation effect of the environment on the 
dynamics of the quantum system which results in the residual coherence of 
the quantum system in the long-time limit. Therefore, it is understandable 
that the non-Markovian effect is a prerequisite for the coherence preservation 
in the long-time limit. Then a natural question is: under the non-Markovian 
dynamics, what is the condition for the quantum coherence to be preserved 
in the long-time limit? This reminds us to examine the eigen solution of the 
whole system, which actually determines the long-time behaviors of quantum 
system. Firstly, let's consider the special case: the environment contains only 
one mode, which corresponds to the J-C model. The whole system possesses 
two real eigenvalues for each excitation-number subspace. Consequently, the 
quantum coherence of the two-level atom would experience loseless oscilla- 
tions when the degree of freedom of the single-mode environment is traced 
out. However, when the environment possesses infinite modes, this would 
not be the truth anymore. It was shown that the real eigenvalues are not 
available anymore (except for the trivial ground state with eigenvalue being 
zero) and the complex eigenvalues are present when the environment has 
infinite modes [35l [38] . This is understandable based on the fact that the 
decay behavior of the quantum system under decoherence is just the effect 
taken by the imaginary part of the complex eigenvalues on the dynamics of 
reduced quantum system. However, John et al found that there is a real 
eigenvalue available when the environment is a PEG medium [35]. Physi- 
cally, the existence of a real eigenvalue means the formation of atom-photon 
bound state. Due to the formation of bound state, spontaneous emission 
would be suppressed and a large proportion of quantum coherence may be 
preserved in the long-time limit. And Ref. [28] reported that the entangle- 
ment can be preserved with a large proportion in the long-time limit when 
two atoms are embedded in the PEG mediums. We argue that the sup- 
pression of the spontaneous emission and the entanglement preservation are 
both contributed from the formation of the bound state. Is this bound state 
available only for such PEG environment or for any environment? Under 
what condition the coherence or entanglement preservation is available for 
generic environments? These questions motivate us to do the investigation 



in Chapters [2] and 131 

There is always lots of entanglement of the quantum system lost irrespec- 
tive of the entanglement could be (partially) preserved or not. Then a nature 
question is: where does the lost entangle go? Modeling the whole systems 
as double J-C model, authors in Ref. |39] have shown that the entanglement 
oscillates between the atoms and the cavities in a lossless way. This is under- 
standable since there is no decoherence in the J-C model. Via introducing a 
normalized collective state, authors in Ref. jlQ] showed that the initial entan- 
glement between the qubits flows entirely to their local environments under 
the Markovian dynamics. We argue that things would be completely different 
under the non-Markovian decoherence dynamics. This judgement is based 
on the following observations. Firstly, it is possible to preserve some entan- 
glement in the quantum system under the non-Markovian dynamics, which 
means that not all of the entanglement between the subsystems is transferred 
to their environments. Secondly, the entanglement preservation is due to the 
formation of the bound state between each subsystem and its environment. 
Therefore, entanglement would exist between each subsystem and its local 
environment. From these facts, we can see that the question where does 
the entanglement go should be reevaluated when the non-Markovian effect is 
taken into account. This motivates us to do the investigation in Chapter |H 

This thesis is organized as follows. In Chapter [2], the decoherence dynam- 
ics of single qubit is studied. We investigate the exact decoherence dynamics 
of a dissipative qubit coupling to a vacuum reservoir. We also study the 
static eigenvalue problem and give the condition when atom-reservoir bound 
state is formed, via which we reveal the mechanism of dynamical decoherence 
suppression due to the bound state. In Chapter [3l entanglement dynamics 
of two qubits under the influence of two independent vacuum reservoirs is 
studied. We give a mechanism of entanglement preservation. In Chapter IH 
we study the entanglement distribution among all possible bipartite parti- 
tions of the same system. Finally, a summary of this thesis and the outlook 
of future works are given in Chapter [51 



Chapter 2 

Decoherence dynamics of a 
dissipative qubit 



In this chapter we study the exact decoherence dynamics of a single qubit 
(two-level atom) in a vacuum reservoir. We compare this result with the one 
obtained under the Markovian approximation. We also study the formation 
of bound state on the decoherence suppression. 

To solve the dynamics for the general open quantum system is rather 
tricky. Here we consider that the environment, with which the qubit inter- 
acts, is in a vacuum state initially. Combining with numerical calculations, 
we can obtain the exact decoherence dynamics of the dissipative qubit. 

2.1 Introduction 

Any realistic quantum system inevitably interacts with its surrounding envi- 
ronment, which leads to the loss of coherence, or decoherence, of the quantum 
system [1]. The decoherence of quantum bit (qubit) is deemed as a main ob- 
stacle to the realization of quantum computation and quantum information 
processing [2] . Understanding and suppressing the decoherence are therefore 
a major issue in quantum information science. For a Markovian environment, 
it is well known that the coherence of a qubit experiences an exponential de- 
crease P^. To beat this unwanted degradation, many controlling strategies, 
passive or active, have been proposed [HI |121 SSI SH US]- 

In recent years much attention has been paid to the non-Markovian effect 



on the decoherence dynamics of open quantum system [16], HTJ HH HHl [501 
[STj [52| [53] . The significance of the non-Markovian dynamics in the study of 
open quantum system is twofold, i) It is of fundamental interest to extend 
the well-developed methods and concepts of Markonian dynamics to non- 
Markovian case [H [12] for the open quantum system in its own right, ii) 
There are many new physical situations in which the Markovian assumption 
usually used is not fulfilled and thus the non-Markovian dynamics has to 
be introduced. In particular, many experimental results have evidenced the 
existence of the non-Markovian effect [TH [191 [201 [21] , which indicates that 
one can now approach the non-Markovian regime via tuning the relevant 
parameters of the system and the reservoir. The non-Markovian effect means 
that the environment, when its state is changed due to the interaction with 
the quantum system, in turn, exerts its dynamical influence back on the 
system. Consequently one can expect decoherence dynamics of the quantum 
system could exhibit a dramatic deviation from the exponential decaying 
behavior. In 2005, DiVincenzo and Loss studied the decoherence dynamics 
of the spin-boson model for the Ohmic heat bath in the weak-coupling limit. 
They used the Born approximation and found that the coherence dynamics 
has a power-law behavior at long-time scale |5l], which greatly prolongs the 
coherence time of the quantum system. Such power-law behavior suggests 
that the non-Markovian effect may play a constructive role in suppressing 
decoherence of the system. Nevertheless, in many cases the finite extension of 
the coherence time of the system is not sufficient for the quantum information 
processing, a question arises whether the coherence of the system can be 
preserved in the long-time limit, even partially. Theoretically, the answer 
is positive if the environment has a nontrivial structure. It has been shown 
that some residual coherence can be preserved in the long-time steady state 
when the environment is a periodic band gap material [25| [26| [27] [28] or leaky 
cavity [231 El] • It is stressed that the residual coherence is due to the confined 
structured environment. A natural question is: Whether the coherence of the 
system can be dynamically preserved or not by the non-Markovian effect if 
the environment has no any special structure, e.g., a vacuum reservoir? 
In this chapter, we study the exact decoherence dynamics of a qubit in- 



teracting with a vacuum reservoir and examine the possibihty of decoherence 
suppression using the non-Markovian effect. The main aim of this chapter is 
to analyze if and how the coherence present in the initial state can be trapped 
with a noticeable fraction in the steady state even when the environment is 
consisted of a vacuum reservoir with trivial structure. We show that the non- 
Markovian effect manifests its action on the qubit not only in the transient 
dynamical process, but also in the asymptotical behavior. Our analysis shows 
that the physical mechanism behind this dynamical suppression to decoher- 
ence is the formation of a bound state between the qubit and the reservoir. 
The no-decaying character of the bound state leads to the inhibition of the 
decoherence and the residual coherence trapped in the steady state. A sim- 
ilar vacuum induced coherence trapping in the continuous variable system 
has been reported in [331 |3l]. Such coherence trapping phenomenon provides 
an alternative way to suppress decoherence. This could be realized by con- 
trolling and modifying the system-reservoir interaction and the properties of 
the reservoir [25] by the recently developed reservoir engineering technique 
[551 1561 [57]. 

2.2 The model and exact decoherence dynam- 
ics of the qubit 

We consider a qubit (two-level atom) which interacts with a vacuum quan- 
tized radiation electromagnetic field. The Hamiltonian of the total system 
reads [1] 

H = wqCt+o-- + ^Wfc4afc + ^(^ffccr+flfc + h.c), (2.1) 

k k 

where Uq is the transition frequency and a± is the raising and lowering op- 
erators of the qubit and a\ and a^, respectively, are the creation and anni- 
hilation operators of the k-th mode with frequency Uk of the radiation field. 
The coupling strength between the qubit and the radiation field is given by 



9k = -y^^, . d, (2.2) 

where e^ and d are unit polarization vector of the radiation field and the 
dipole moment of the qubit respectively. Thoughout this paper we assume 



h = i. 

To obtain the exact djTiainics of the qubit, we first consider the following 
two simple cases. For simplicity, we assume there is no correlation between 
the qubit and its reservoir at the initial time t = 0. If the initial state of 
the system is |\E'(0)) = |— ,{0^}), where |— ) denotes the ground state of 
the qubit and 1(0,^}) represents the vacuum state of the reservoir, the whole 
system will not evolve with time under the Hamiltonian f l2.ip . Whereas to 
the initial state |\E'(0)) = |+, {0^}), in which |+) denotes the exited state of 
the qubit, the time evolution of the total system has the following form 

|vl>(t)) = 6o(t) |+, {Ok}) + J2 ^fc W I-' il'^}) ' (2-3) 

k 

where |{lfc}) represents the field state containing one photon in the k-th 
mode. Applying the Schrodinger equation, we get the time evolution of the 
probability amplitudes 

iboit) = boitPo + J^aMt), (2.4) 

k 

ihk{t) = bkit)uk + g*Mt), (2.5) 

where the superscript dot represents the differential with respect to time. 
Solving Eq. f l2.5p formally and substituting the solution into Eq. f l2.4p . we 
can obtain 

kit) + ttuoboit) = - I f{t- r)bo{r)dT, (2.6) 



^0 

where the kernel function is f{x) = YlT=o \9k\ exp(— zcjfco;). Obviously the 
memory effect has been registered in the kernel function. In the continuous 
limit of the environment frequency, the kernel function has the form 

f{x) = / J(w)e-*'^^du;, (2.7) 

Jo 

where J{u}) = rjuj^e'^ is the spectral density [58], which characterizes the 
coupling strength of the reservoir to the qubit with respect to the reservoir 
frequency and r] = ,2wc)^2 — ■ '^° eliminate infinity in frequency integration, 
we have introduced the cutoff frequency Uc- On physical grounds, the intro- 
ducing of the cutoff frequency means that not all of the infinite modes of the 



10 



reservoir contribute to the interaction with the qubit, and one always expects 
the spectral density going to zero for the modes with frequencies higher than 
certain characteristic frequency. It is just this characteristic frequency which 
determines the specific behavior and the properties of the reservoir. One can 
see that in our model, the spectral density has a super-Ohmic form |58] . 

From the time evolution of the above two situations, one can get the time 
evolution of any given initial state of the system readily. For an initially 
mixed state, which is described by the following density operator 

Ptot(O) = (pil| + )(+|+Pl2|+)(-|+P2l|-)(+| 

+P22|-)(-|)®|{0U)({0}fc|. (2.8) 

The time evolution of the total system can be calculated explicitly. In fact, 
what we care about is the reduced density matrix of the qubit, which is 
obtained by tracing over the reservoir variables 

^(+\ - f Pn\Ht)f Pi2bo{t) \ . . 

'^'^-[ P.Mt) l-pnIWrJ- ^'-'^ 

Differentiating Eq. (12.91) with respect to time, we may obtain the equation 
of motion of the qubit 

pit) = -^^Ka_,p(t)] + ^[2a_p(t)a+ 

-a^a^pit) - pit)a+a.], (2.10) 

where fi(t) = -2Im[^] and 7(t) = -2Re[^]. n{t) plays the role of time- 
dependent shifted frequency and •yit) that of time-dependent decay rate [1]. 
It is worth mentioning that during the derivation of master equation (I2.10p 
we have not resorted to the Born-Markovian approximation. Therefore Eq. 
(I2.10p is the exact master equation of the qubit system. 

It is interesting to notice that one can reproduce the conventional Marko- 
vian one from our exact non-Markovian master equation under certain ap- 
proximations. By redefining the probability amplitude as 6o(t) = bQ{T)e~^'^°'^ , 
one can recast Eq. (12.61) into 

^OO ft 

h'^{t)+ dcoJico) dTe'^'^°-''^^'-^^b'f^{T)=0, (2.11) 

Jo Jo 
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where J{oj) is defined the same as above. Then, we take the Markovian 
approximation, 

h',{T)=h',{t). (2.12) 

namely, approximately taking the dynamical variable to the one that depend 
only on the present time so that any memory regarding the earlier time 
is ignored. The Markovian approximation is mainly based on the physical 
assumption that the correlation time of the reservoir is very small compared 
with the typical time scale of system evolution. Also under this assumption 
we can extend the upper limit of the r integration in Eqs. (12.111) to infinity 
and use the equality 

lim f cire^^('^o-a;)(t-r) _ ^^(^ _ ^^^ ^ -pf — \ — y ^2.13) 

t-s>oo Jq \uJ — UJq/ 

where P and the delta-function denote the Cauchy principal value and the 
singularity, respectively. The integro-differential equation in (12. lip is thus 
reduced to a linear ordinary differential equation. The solutions of 6q as well 
as bo can then be easily obtained as 

where 5co = P [^ ^^' ^ . Thus one can verify that, 

7(t) = 70 = 27rJ(a;o), ^(t) = ^o = 2(a;o - f^w), (2.15) 

which are exactly the coefficients in the Markovian master equation of the 
two-level atom system p. 

2.3 Purity and decoherence factor 

To quantify the decoherence dynamics of the qubit, we introduce the follow- 
ing two quantities. The first one is the purity, which is defined as [2] 

Vit) = Trp2(t). (2.16) 

Clearly p = 1 for pure state and p < 1 for mixed state. The second quantity 
describing the decoherence is the decoherence factor c(t) of the qubit, which 
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is determined by the off-diagonal elements of the reduced density matrix 

|pi2(t)|=c(t)|pi2(0)|. (2.17) 

The decoherence factor maintains unity when the reservoir is absent and 
vanishes for the case of completely decoherence. 

For definiteness, we consider the following initial pure state of the qubit 

|^(0)) = a|+)+/3|-), (2.18) 

in which a and /3 satisfy the normalization condition. Using Eq. (12. 9p . the 
exact time evolution of the qubit is easily obtained 

m = ( I"': I'ff , f ??f> , ) . (2.19) 

V ctfihl{t) 1- |a| %{t)\ ) 

With Eq. (12.191) . the purity and decoherence factor can be expressed explic- 
itly 

V{t) = 2\at%{tt{%(t)f -\\^\ (2.20) 

and 

c(t) = |6o(t)|. (2.21) 

It is easy to verify, under the Born-Markovian approximation, the purity 
and decoherence factor have the following forms 

p(t) = 1\a\^ e-^«*(e-^o* - 1) + 1, (2.22) 

and 

c(t) = e-^*, (2.23) 

where the time-independent decay rate 70 is given in Eq. (12.151) . Obvi- 
ously, the system asymptotically loses its quantum coherence (c(oo) = 0) 
and approaches a pure steady state {p{oo) = 1) irrespective of the form of 
the initial state under the Markovian approximation. One can also find from 
Eqs. (12. 20112. 231) that the probability amplitude of excited state plays key 
role in the decoherence dynamics. 
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Figure 2.1: Time evolution of 7(t), p{t) and c{t) in non-Markovian situation 
(solid line) and the corresponding Markovian situation (dashed line), when 
r/ and Uc are small. The parameters used here are a = 1/a/2, t] = 0.08 and 
cjjujo = 1.0. 

2.4 Numerical results and analysis 

In this section, by numerically solving Eq. (12. 6p . we study the influence of 
memory effect of reservoir on the exact dynamics of the qubit. Noticing the 
fact that the memory effect registered in the kernel function is essentially 
determined by the spectrum density J^u), one can expect that J{ou) plays 
an major role in the exact dynamics of the qubit. In the following, we show 
how the decoherence of the qubit can be fully suppressed under the non- 
Markovian dynamics in terms of the relevant parameters of J{uj) [59] . 

2.4.1 The influence of coupling constant 

In the following, we numerically analyze the exact decoherence dynamics of 
the qubit with respect to decay rate •yit), purity p(t) and decoherence factor 
c{t) in terms of the coupling constant r] [59] . 

In Fig. 12.11 we plot the time evolution of decay rate 7(t), purity p{t), 
decoherence factor c(t) and their Markovian correspondences in the weak 
coupling and low cutoff frequency case. We can see that 7(t) shows distinct 
difference from its Markovian counterpart over a very short time interval. 
With time, 7(t) tends to a deflnite positive value. The small "jolt" of 7(t) 
in the short time interval just evidences the backaction of the memory effect 
of the reservoir exerted on the qubit [GOj. It manifests that the reservoir 
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Figure 2.2: Time evolution of 7(t), p{t) and c{t) in non-Markovian situation 
(solid line) and the corresponding Markovian situation (dashed line), when r) 
is large. The parameters used here are a = l/-\/2, t] = 1.0 and cOc/coq = 1.0. 

does not exert decoherence on the qubit abruptly, just as the result based on 
Markovian approximation, but dynamically influences the qubit and grad- 
ually establishes a stable decay rate to the qubit. Furthermore, it is also 
shown that the decay rate is positive in the full range of evolution, which 
results in any initial qubit state evolving to the ground state \i/j{oo)) = |— ) 
irreversibly. Consequently the decoherence factor monotonously decreases to 
zero with time and the purity approaches unity in the long-time limit, which 
is consistent with the result under Markovian approximation. The result 
indicates that although the reservoir has backaction effect on the qubit, it 
is quite small. And the dissipation effect of the reservoir dominates the dy- 
namics of the qubit. Thus no qualitative difference can be expected between 
the exact result and the Markovian one with the backaction effect ignored. 
Therefore the widely used Markovian approximation is applicable in this case. 
Nevertheless, at the short and immediate time scales the overall behavior is 
still quite different from that of the Markovian dynamics. The decoherence 
factor shown in the righ-hand panel of Fig. 12. II shows non-exponential decay, 
which is in agreement with the result obtained previously in the spin-boson 
model in the weak-coupling limit [M]- However, the situation is dramatically 
changed if the coupling is strengthened as discussed below. 

With the same cutoff frequency as in Fig. 12.11 but a larger coupling 
constant, we plot in Fig. 12.21 the decay rate, purity and decoherence factor 
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in the strong coupling case. In this case the non-neghgible backaction of 
the reservoir has a great impact on the dynamics of the qubit. Firstly, we 
can see that the decay rate not only exhibits oscillations, but also takes 
negative values in the short time scale. Physically, the negative decay rate 
is a sign of strong backaction induced by the non-Markovian memory effect 
of the reservoir. And the oscillations of the decay rate between negative 
and positive values reflect the exchange of excitation back and forth between 
qubit and the reservoir [Hj. Consequently both the decoherence factor and 
the purity exhibit oscillations in a short-time scale, which shows dramatic 
deviation to the Markovian result. Therefore, entirely different to the weak 



coupling case in Fig. 12. ![ the reservoir in the strong coupling case here has 
strong backaction effect on the qubit. Secondly, we also notice that the decay 
rate approaches zero in the long-time limit. The vanishing decay rate means, 
after several rounds of oscillation, the qubit ceases decaying asymptotically. 
The non-Markovian purity maintains a steady value asymptotically, which is 
less then unity. This indicates that the steady state of the qubit is not the 
ground state anymore, but a mixed state. The decoherence factor also tends 
to a non-zero value, which implies that the coherence of the qubit is preserved 
with a noticeable fraction in the long-time steady state. These phenomena, 
which are qualitatively different to the Markovian situation, manifest that 
the memory effect has a considerable contribution not only to the short-time, 
but also to the long-time behavior of the decoherence dynamics. The presence 
of the residual coherence in the steady state also suggests a potential active 
control way to protect quantum coherence of the qubit from decoherence via 
the non-Markovian effect. 

2.4.2 The influence of cutoff frequency 

The cutoff frequency Uc, on the one hand, is introduced to eliminate the 
infinity in the frequency integration. On the other hand it also determines 
the frequency range in which the power form is valid [61]. In the following, we 
elucidate the influence of cutoff frequency on the exact decoherence dynamics 
[59]. 

Fixing 1] as the value in Fig. 12.11 and increasing the cutoff frequency, we 
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Figure 2.3: Time evolution of 7(t), p{t) and c{t) in non-Markovian situation 
(solid line) and the corresponding Markovian situation (dashed line), when c^c 
is large. The parameters used here are a = 1/a/2, i] = 0.08 and cOc/coq = 3.0. 

plot in Fig. 12.31 the dynamics of the qubit in a high cutoff frequency case. 
It shows that a similar decoherence behavior as the strong coupling case in 
Fig. 12.21 can be obtained. After several rounds of oscillation, the decay rate 
tends to zero in the long-time limit. The negative decay rate makes the 
lost coherence partially recovered. The vanishing decay rate in the long-time 
limit results in the decoherence frozen before the qubit gets to its ground 
state. Thus there is some residual coherence trapped in the steady state. 
Similar to the strong coupling case, it is essentially the interplay between 
the backaction and the dissipation on the dynamics of qubit which results in 
the inhibition of decoherence. We argue that in this high cutoff frequency 
regime, the widely used Markovian approximation is not applicable because 
of the strong backaction effect of the reservoir. 

2.4.3 The physical mechanism of the decoherence in- 
hibition: the formation of atom-photon bound 
state 

From the analysis above we can see clearly that the decoherence can be 
inhibited in the non-Markovian dynamics. A natural question is: What is 
physical mechanism to cause such dynamical decoherence inhibition? To 
answer this question, let us find the eigen solution of Eq. (12. ip in the sector 
of one-excitation in which we are interested ^9\. The eigenequation reads 
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H\'^e) = E\(fE), where \ipe) = Co|+,{Ofc}) + ZlfcLo^fc h, Ifc)- After some 
algebraic calculation, we can obtain a transcendental equation of E 

y{E) ^uo- r ^^du = E. (2.24) 

Jo tJ - E 

From the fact that y{E) decreases monotonically with the increase of E when 
E < we can say that if the condition y{0) < 0, i.e. 

uo - 2r]^ < (2.25) 

is satisfied, y{E) always has one and only one intersection in the regime E < 
with the function on the right-hand side of Eq. (I2.24p . Then the system will 
have an eigenstate with real (negative) eigenvalue, which is a bound state 
[551 EH], in the Hilbert space of the qubit plus its reservoir. While in the 
regime of E' > 0, one can see that y{E) is divergent, which means that no 
real root E can make Eq. (I2.24p well-defined. Consequently Eq. (12.241) 
does not have positive real root to support the existence of a further bound 
state. It is noted that Eq. (I2.24p may possess complex root. Physically 
this means that the corresponding eigenstate experiences decay contributed 
from the imaginary part of the eigenvalue during the time evolution, which 
causes the excited-state population approaching zero asymptotically and the 
decoherence of the reduced qubit system. 

The formation of bound state is just the physical mechanism responsible 
for the inhibition of decoherence. This is because a bound state is actually a 
stationary state with a vanishing decay rate during the time evolution. Thus 
the population probability of the atomic excited state in bound state is con- 
stant in time, which is named as "population trapping" [25l[27]. This claim 
is fully verified by our numerical results. The parameters in Fig. 12.11 do not 
satisfy the condition (12.251) to support the existence of a bound state, then 
the dynamics experiences a severe decoherence. While with the increase of 
either rj (in Fig. 12. 2p or Uc (in Fig. 12.30 . the bound state is formed. Then 
the system and its environment is so correlated that it causes the decay 
rate of the system in the non-Markovian dynamics exhibiting: 1) transient 
negative value due to the backaction of the environment; 2) vanishing asymp- 
totic value. Such interesting phenomenon, i.e. the vanishing asymptotical 
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Figure 2.4: Time evolution of c(t) in the non-Markovian dynamics with dif- 
ferent 7] when Uc/ojq = 1.0 (upper panel) and with different Uc when rj = 0.08 
(lower panel). 



decay rate in the large cutoff frequency regime for super-Ohmic spectrum 
density, was also revealed in Ref. |62]. This effect of course is missing in 
the conventional Born-Markovian decoherence theory, where the reservoir is 
memory less. 

In order to understand the exact decoherence dynamics more completely, 
we plot in Fig. 12.41 the crossover from coherence destroying to coherence 
trapping via increasing either the coupling constant or the cutoff frequency. 
Coherence trapping can be achieved as long as the bound state is formed. 
Therefore, one can preserve coherence via tuning the relevant parameters of 
system and the reservoir, e.g. the qubit-reservoir coupling constant and the 
property of the reservoir so that the condition (12.251) is satisfied. 
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2.5 Summary 

In summary, we have investigated the exact decoherence dynamics of a qubit 
in a dissipative vacuum reservoir. We have found that even in a vacuum en- 
vironment without any nontrivial structure, we can still get the decoherence 
suppression of the qubit owing to the dynamical mechanism of the non- 
Markovian effect. From our analytic and numerical results, we find that the 
non-Mar kovian reservoir has dual effects on the qubit: dissipation and back- 
action. The dissipation effect exhausts the coherence of the qubit, whereas 
the backaction one revives it. In the strong coupling and/or high cutoff fre- 
quency regimes, a bound state between the qubit and its reservoir is formed. 
It induces a strong backaction effect in the dynamics because the reservoir is 
strongly correlated with the qubit in the bound state. Furthermore, because 
of the non-decay character of the bound state the decay rate in this situation 
approach zero asymptotically. The vanishing of the decay rate causes the 
decoherence to cease before the qubit decays to its ground state. Thus the 
qubit in the non-Markovian dynamics would evolve to a non-ground steady 
state and there is some residual coherence preserved in the long-time limit. 
Our results make it clear how the non-Markovian effect shows its effects on 
the decoherence dynamics in different parameter regimes. 

The presence of such coherence trapping phenomenon actually gives us 
an active way to suppress decoherence via non-Markovian effect. This could 
be achieved by modifying the properties of the reservoir to approach the 
non-Markovian regime via the potential usage of the reservoir engineering 
technique [551 EH [571 [23]. Many experimental platforms, e.g. mesoscopic 
ion trap [55|, [56], cold atom BEG [57], and the photonic crystal material 
[25] have exhibited the controllability of decoherence behavior of relevant 
quantum system via well designing the size (i.e. modifying the spectrum) 
of the reservoir and/or the coupling strength between the system and the 
reservoir. It is also worth mentioning that a proposal aimed at simulating 
the spin-Boson model, which is relevant to the one considered in this paper, 
has been reported in the trapped ion system [6^. On the other side many 
practical systems can now be engineered to show the novel non-Markovian 
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effect [m [201 EB 123 • ^11 these achievements show that the recent advances 
have paved the way to experimentally simulate the paradigmatic models of 
open quantum system, which is one part of the new-emergent field, quantum 
simulators [65] . Our work sheds new light on the way to indirectly control and 
manipulate the dynamics of quantum system in this experimental platforms. 
A final remark is that our results can be generalized to the system con- 
sisted of two qubits, each of which interacts with a local reservoir. Because 
of the coherence trapping we expect that the non-Markovian effect plays 
constructive role in the entanglement preservation [281 EH EH] . 
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Chapter 3 

Mechanism of entanglement 
preservation 



In this chapter, we study the mechanism of entanglement preservation. We 
found entanglement can be preserved in the long time limit as long as bound 
states are formed in the local systems. We also find the non-Markovian effect 
has a profound effect on the entanglement preservation. 

3.1 Introduction 

Entanglement is not only of fundamental interest to quantum mechanics, but 
also of great importance to quantum information processing [2]. However, 
due to the inevitable interaction of qubits with their environments, entangle- 
ment always experiences degradation. A phenomenon that the entanglement 
between two qubits may completely disappear at a finite time, known as "en- 
tanglement sudden death" (ESD), has been predicted theoretically [121 [S] 
and subsequently been verified experimentally [HI [17] , which indicates the 
specific behavior of the entanglement different from the coherence. From the 
point of view of applications, the ESD is apparently disadvantageous to the 
quantum information processing. 

Recently, Bellomo et al. [29j found that the entanglement can revive after 
some time interval of the ESD and thus extends significantly the entangled 
time of the qubits. This remarkable phenomenon, which has been experimen- 
tally observed [22], is physically due to the dynamical backaction, i.e., the 
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non-Markovian effect, of the memory environments. However, in many cases 
the finite extension of the entangled time is not enough and thus it is desired 
to preserve a significant fraction of the entanglement in the long time limit. 
Indeed, it was shown [28] that some noticeable fraction of entanglement can 
be obtained by engineering structured environment such as photonic band- 
gap materials [35]. According to these works, it is still unclear if the residual 
entanglement is fundamentally due to the specific structured materials or due 
to certain physical mechanism. Is there any essential relationship between 
the ESD and/or its revival phenomena and the residual entanglement? 

In this chapter we focus on these questions and elucidate the physical 
nature of the residual entanglement. Before proceeding, it is helpful to recall 
the physics of quantum electrodynamics of a single two-level atom placed in 
a dielectric with a photonic band gap [26| [35] . The coupling between the 
excited atom and electromagnetic vacuum in the dielectric leads to a novel 
photon-atom bound state, in which the fractional atomic population on the 
excited state occurs, also known as population trapping [27]. This result has 
been verified experimentally for quantum dots embedded in photonic band- 
gap environment [37]. The population trapping has been directly connected 
to the entanglement trapping due to the structured environment [28]. Here 
we reveal for the first time that there are two essential conditions to pre- 
serve the entanglement. One is the existence of the bound state between the 
system and its environment, which provides an ability to preserve the entan- 
glement, and the other is the non-Markovian effect, which provides a way 
to preserve the entanglement. Our result can reproduce the ESD [13] when 
the non-Markovian effect is neglected. The phenomenon of the ESD and its 
revival discussed in Ref. [29] results from the non-Markovian effect when the 
bound state is not available. The interplay between the availability of the 
bound state and the non-Markovian effect can lead to a significant fraction 
of the entanglement preserved in the steady state. We verify these results 
by considering two reservoirs modeled by the super-Ohmic and Lorentzian 
spectra, respectively. The result provides a general method on how to protect 
the entanglement by engineering the environment. 
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3.2 The model and entanglement dynamics 

We now consider two spatially separated systems A and B, each has a two- 
level atom coupled to a vacuum reservoir, and the two qubits are initially 
entangled but have no direct interaction. Owing to the independence of the 
two systems |29J, we can investigate single "qubit + reservior" system at the 
first place, then extend our studies to the double-one. 

The single "qubit + reservoir" system can be formulated by the following 
Hamiltonian 

H = u]oa+a_ + ^ cufe&l&fc + ^igk(y+h + glf^-bl), (3.1) 

k k 

where ujq is the transition frequency of the two-level atom, and a± are the 
atom raising and lowering operators, h\ and hk are respectively the creation 
and annihilation operators of the k-th mode with frequency Uk of the reser- 
voir. Qk denotes the coupling strength between the atom and the radiation 
field. 

Following the procedure we done in the Section 2.2, we can obtain the 
master equation of the qubit, 

dp^it) .Vtit)^ ^, ^^ 



dt 2 

■[2a_p^(t)a+ - cj+a_p'{t) - p^'(i)a+(T_], (3.2) 



7(^)ro, ^su\^ ^ ^ ^S( 



where Vt{t) = — 2Im[^^], and 7(t) = — 2Re[^44]. This is exactly the master 
equation of the single qubit. fl{t) and 7(t) play, respectively, the role of 
time-dependent Lamb shift and decay rate. 

With the dynamics of single "qubit + reservoir " , we can readily study 
the decoherence dynamics of the double-one. We assume, for simplicity, the 
two systems are the same, and the double-system is initially in a mixed state 



pL(o) 



f Pll Pl2 Pl3 Pl4 ^ 

Pl2 P22 P23 P24 

Pl3 P23 P33 P34 

\ Pl4 P24 P34 P44 / 



®l{0}n)({0}n|. (3.3) 



24 



Following the method we applied in the single qubit system, it is easy to 
obtain the time evolution of the two qubits. The diagonal elements are 

Piiit) = Pu\co{t)\^, 

P22it) = P22|Co(t)|Vpil|Co(t)|^(l- |Co(t)|^), 
Pisit) = P33|Co(t)|^+Pll|Co(t)|^(l- |Co(t)|^), 

pUt) = l + Pii|co(t)|^-|co(t)|'(2pii+p22 + P33), (3.4) 

the nondiagonal elements are 

Puit) = Pi2\co{t)fco{t), 
Pnit) = Pi3\coit)fco{t), 

PuW = Pl4Co(t), 

pLit) = P23\co{t)f, 

P24W = P24Co(t) +Pl3Co(i)(l - |co(t)|^), 

Puit) = P34Co(t)+Pl2Co(t)(l-|co(t)|'), (3.5) 

and pfj(t) = pj*(t). Differentiating the reduced density matrix with respect 
to time, we can obtain the equation of motion for the two qubits 



dt 2 

.lit) 



.^([a>f,p^(t)] + [a«a«,p^(t)] 

{[2a^p^{t)al - a^a^p^it) - p^(t)a^a^] 



2 
+ [2a^p^it)a^ - a^a^_p^{t) - p^{t)a^a^]}, (3.6) 

where Q{t) and 7(t) are defined the same as before. 

To investigate the entanglement dynamics of the bipartite system, we 
apply Wootters concurrence [67]. The concurrence can be calculated ex- 
plicitly from the time dependent density matrix p'^{t) of the two qubits, 
C{p^) = max{0, a/Ai — V^ — y/Xs — y/X^}, where the quantities Aj are the 
eigenvalues of the matrix C 

C = pVX)p^*«®0, (3.7) 

arranged in decreasing order. Here p'^* means the complex conjugation of 
P"^, and ay is the Pauli matrix. It can be proved that the concurrence varies 
from for a separable state to 1 for a maximally entangled state. 
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As is pointed out, in the Markovian situation the ESD occurs due to the 
double excitation of the initial state in a vacuum reservoir [15]. In what 
follows, we consider the initial state of which the concurrence dynamics can 
exhibit ESD in Markovian approximation, and compare this with the non- 
Mar kovian situation. 

For an initially entangled pure state in the standard bases 

^(0) = a| )+/3|++), (3.8) 

a and /3 satisfy normalization condition. From the time-dependent reduced 
density matrix of the two qubits and the definition of concurrence, we obtain 

C(p^) = max{0, 2 \c^{tt |/3| [|«| - |/3| (1 - |co(t)|')]}. (3.9) 

From the expression of concurrence, it can be found that the behavior 
of time-dependent factor local excited state population (|co(t)| ) completely 
determines the dynamics of concurrence. In particular, if local system deco- 
herences completely (|co(oo)| = 0), entanglement would vanish. Whereas, if 
decoherence in local system is inhibited, then it is possible to protect entan- 
glement in the long-time limit. 

3.3 Mechanism of entanglement preservation 

Entanglement is the fundamental resource of quantum information process- 
ing. Entanglement preservation is a key step for entanglement applications. 
To achieve this aim, as we have discussed in Sec. 3.2, one needs the de- 
coherence suppression in the local system. This could be obtained when 
atom-photon bound state is formed. 

Following the same procedure done in Sec. 2.4.3, we obtain the condition 
for the formation of bound state [59l [66] 



y{E) ^000- r ^^du = E, (3.10) 

Jo tJ - E 

where J(w) = J2k lfi'fcP^('^ '~ ^k) is the spectral density of the reservoir. The 
solution of Eq. (13.101) highly depends on the explicit form of the J{(jj). If 
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the reservoir contains only one mode w', then J{uj) = g'^S{uj — u'). This is 
the ideal Jaynes-Cummings model [68], in which two bound states in the 
one excitation sector are formed and as a result the dynamics of the system 
displays a lossless oscillation. When the reservoir contains infinite modes, one 
can model J{u}) by some typical spectrum functions such as the super-Ohmic 
or Lorentzian form. 

We firstly consider the super-Ohmic spectrum J(uj) = 'q^e^'^^^", where 
?7 is a dimensionless coupling constant and Uc characterizes the frequency 
regime in which the power law is valid [61j. It corresponds to that the 
reservoir consists of a vacuum radiation field, where gk oc ^yuJk j^. The 
existence of a bound state requires that Eq. fl3.1Up has at least a real solution 
for E < 0. It is easy to check that the solution always exists if the condition 
y{0) < 0, i.e. Uq — 2ri^ < is satisfied. Otherwise, no bound state exists. 
This condition can be fulfilled easily by engineering the environment. For the 
Lorentzian spectrum it is found that a criterion when a bound state exists 
can not be obtained analytically. In this case one can use the diagrammatic 
technique shown later. 

The formation of bound state may lead to the inhibition of spontaneous 
emission and results in the population trapping. Therefore the formation of 
local bound state has a profound effect on the preservation of non-local co- 
herence(entanglement). One could imagine that non-Markovian effect is also 
a key factor to preserve entanglement since non-Markovian effect is a kind 
of backaction effect which could compensate the lost coherence of quantum 
system. In the following, we explicitly study the formation of bound state 
and non-Markovian effect on entanglement dynamics using the examples of 
super-ohmic and Lorentzian vacuum reservoir |66j . 

Consider firstly the super-Ohmic case. Fig. 13.11 shows the entanglement 
dynamics in different parameter regimes, i.e., {uc, rf) = (O.Tcjq, 0.2), (O.Twq, 1.0) 
and (S.Owo, 0.2). For the first two parameter sets the bound state is absent, 
while for the last one it is available, as shown in Fig. 13.11 (a). Whether the 
bound state exists or not plays a key role in the entanglement preservation 
in the long time limit. When the bound state is absent, the residual entan- 
glement approaches zero in a long enough time, as shown by the solid lines in 
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Figure 3.1: Entanglement dynamics of the two-qubit system with local 
super-Ohmic reservoirs, (a) Diagrammatic solutions of Eq. f l3.10p with dif- 
ferent parameters. C{t) as a function of time is shown in (b): {uJc,f]) = 
(0.7^0,0.2), (c): (ucv) = (0.7wo, 1-0) and (d): {uc,r]) = (B-Ocuq, 0.2). The 
parameter a is taken as 0.7. For comparison, C{t) under the Markovian 
approximation has also been presented by using the same parameters. 
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Figure 3.2: The decay rate r(t) as a function of time in the non-Markovian 
and Markovian cases. The parameters used are Uc = S.Oojq and r] = 0.2. 



Fig. 13.11 (b) and (c). Difference between these two cases is that Fig. 13.1( b) 
is in weak couphng regime, where the non-Markovian effect is weak, while 
Fig. 13.1( c) is in strong coupling regime, where the strong non-Markovian 
effect leads to an obvious oscillation. When the bound state is available, 
the situation is quite different, as shown in Fig. 13.1( d). The entanglement 
firstly experiences some oscillations due to the energy and/or information 
exchanging back and forth between the qubit and its memory environment 
[IH], then approaches a definite value in the long time limit, where the decay 
rate approaches zero after some oscillations, as shown in Fig. 13.21 The entan- 
glement preservation is a result from the interplay between the existence of 
the bound state (providing an ability to preserve the entanglement) and the 
non-Markovian effect (providing a way to preserve the entanglement) |66j . 
The claim can be further verified by that the entanglement preservation is 
absent in the Markovian dynamics, as shown by the dashed lines in Fig. 13.11 
(b), (c) and (d), where the entanglement displays sudden death irrespective 
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Figure 3.3: The residual entanglement for different initial states with a = 
0.7,0.5 and 0.3. The other parameters used are the same as those in Fig 
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of the availability of the bound state. It is because the Markovian environ- 
ment has no memory and the energy/information flowing from the qubit to 
its environment is irreversible and the decay rate keeps to be fixed (see, Fig. 
13.21). In this case one has 



C{t) = max{0, 2e~'^»* \/3\ [\a\ - \/3\ (1 



-2rnt 



n]}, 



(3.11) 



which shows a finite disentanglement time when |a| < |/3| [14j. In a word, the 
above discussion manifests clearly two conditions to preserve the entangle- 
ment, i.e., the availability of the bound state and the non-Markovian effect 
|66j . not only the structured environment as emphasized in Ref. [28] . 

The above discussion focused on almost maximally entangled initial state 
by taking a = 0.7. In Fig. 13.31 we show the results for different initial states 
with different initial entanglement. With decreasing the initial entanglement, 
the residual entanglement also decreases in the long time limit and finally, 
the ESD happens for a = 0.3. The result can be understood from Eq. (13. 9p . 
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Figure 3.4: The entanglement dynamics with the Lorentzian spectrum, 
(a) Diagrammatic solutions of Eq. fl3.10p with different parameters A = 
O.lwo, 2.0a;o and ISwo- (b) C(t) as a function of time for the corresponding 
three parameter regimes. The insert in (b) are the decay rate as a function 
of time. The other parameters used are 7 = S.Owq and a = 0.7. 

On the one hand, the residual entanglement is determined by co(oo), which 
is directly related to the property of the bound state. On the other hand, 
the residual entanglement is also determined by the competition between 
the first and the second terms in Eq. (13. 9p , which is dependent of the initial 
state. 

In order to make a comparative study and confirm our observations we 
consider the Lorentzian spectrum if the reservoir is composed of lossy cavity. 



Au) 



1 



7A2 



2n (w - uo)^ + A2 



(3.12) 



where 7 is the coupling constant and A is the spectrum width. This model 
has also been studied in Ref. [29], where the lower limit of frequency integral 
in f{t — r) was extended from zero to negative infinity. This extension is 
mathematically convenient but the availability of the bound state is missed. 
Here we follow the original definition of the frequency integral ranges. 

Our model with the Lorentzian spectral density corresponds exactly to 
the extended damping J-C model [I]. It is noted that the strong coupling 
of J-C model has been achieved in circuit QED |69] and quantum dots [70] 
systems. Fig. 13.41 shows the entanglement dynamics of the qubits under 
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the Lorentizan reservoir for different spectrum widths in the strong couphng 
regime. When A = O.lwo, Eq. fl3.10p lacks the bound state. According to the 
above discussion, there is no residual entanglement in the long time limit. 
This is indeed true, as shown in Fig. 13.41 (b). However, it is noted that 
before becoming zero the entanglement exhibits "sudden death" and revives 
after some time for several times. This is an analog of the central result 
found in Ref. [2H], i-e., the phenomenon of the ESD and revival. Apparently, 
this is due to the non-Markovian effect, the revival is a result of backaction 
of the memory reservoir. When increasing the spectrum width the bound 
states become available, the situation changes. The significant fraction of 
the entanglement initially present is preserved in the long time limit, where 
the decay rates shown in the insert of Fig. 13.41 (b) approach zero in these 
cases. Likewise, the physical nature of the entanglement preservation is still 
the interplay between the bound state and the non-Markovian effect. The 
more stronger the coupling is, the more striking the entanglement oscillates 
as a function of time, consequently, the more noticeable the non-Markovian 
effect is, as shown in Fig. 13.51 For 7 = 0.2a;o5 the system is in the weak 
coupling regime, where the bound state is also not available. As a result, the 
ESD is reproduced in this case. 
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In Fig. I3.6l we present a phase diagram of the entanglement in the steady 
state for the Lorentzian spectral density. In the large 7 and small A regime. 



32 







With Residual Entanglement 



Without Residual Entanglement 







5 10 

A,/co„ 



15 



Figure 3.6: The phase diagram of the residual entanglement in the steady 
state for Lorentzian spectrum. 



the system approaches the J-C model. In this situation the strong backaction 
effect of the reservoirs makes the qubit system hard to form a steady state. 
The entanglement oscillates with time but has no dissipation. In the small 7 
and large A regime, the non-Markovian effect is extremely weak and our re- 
sults reduce to the Markovian case. In a limit of the flat spectral density, the 
Born-Markovian approximation is applicable and the system has no bound 
state. This is the case of the ESD HH [661. 



3.4 Summary 

In this chapter, we have studied the entanglement protection of two qubits 
in two uncorrelated reservoirs. Two essential conditions to preserve the en- 
tanglement are explored, one is the existence of the bound state of the sys- 
tem and its reservoir and the other one is the non-Markovian effect. The 
bound state provides the ability of the entanglement preservation and the 
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non-Markovian effect provides the way to protect the entanglement. The 
previous results on the entanglement dynamics in the literature can be con- 
sidered as the specific cases where these two conditions have not been fulfilled 
at the same time. The result provides a unified picture for the entanglement 
dynamics and gives a clear way how to protect the entanglement. This is 
quite significant in the quantum information processing. 

The presence of such entanglement preservation gives us an active way to 
suppress decoherence. This could be achieved by modifying the spectrum of 
the reservoirs to approach the non-Markovian regime and form a bound state 
via the potential usage of the reservoir engineering technique [55l [TH 172] . 
Fortunately, we notice that many practical systems have now be engineered 
to show strong non-Markovian effect fi9[ [20t [211 [22]. AH these achievements 
have paved the way to experimentally simulate the paradigmatic models of 
open quantum system, which gives a hopeful prospective to preserve the 
entanglement. 
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Chapter 4 

Entanglement distribution and 
its invariance 



In this chapter we study the entanglement distribution among bipartite sys- 
tems of quantum systems and their reservoirs. Following the method used in 
Ref. [lO], we find an invariant and entanglement can be distributed among 
all the bipartite subsystems. 

4.1 Introduction 

In Chapter 3, we have study the entanglement dynamics under different envi- 
ronments. We have figured out the mechanism of entanglement preservation, 
i.e. entanglement can be preserved when bound states are formed under the 
non-Markovian dynamics. 

Recently, Lopez et al. asked a question about where the lost entanglement 
between the qubits goes |1D]. Interestingly, they found that the lost entan- 
glement of the qubits is exclusively transferred to the reservoirs under the 
Markovian dynamics and the ESD of the qubits is always accompanied with 
the entanglement sudden birth (ESB) of the reservoirs. This means that the 
entanglement does not go away, it is still there but just changes the location. 
This is reminiscent of the work of Yonac et al. [39] , in which the entangle- 
ment dynamics has been studied in a double J-C model. They found that the 
entanglement is transferred periodically among all the bipartite partitions of 
the whole system but an identity (see below) has been satisfied at any time. 
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This may be not surprising since the double J-C model has no decoherence 
and any initial information can be preserved in the time evolution. However, 
it would be surprising if the identity is still valid in the presence of the de- 
coherence, in which a non-equilibrium relaxation process is involved. In this 
chapter, we show that it is indeed true for such a system consisted of two 
qubits locally interacting with two reservoirs. We find that the distribution 
of the entanglement among the bipartite subsystems is dependent of the ex- 
plicit property of the environment and its coupling with the qubit. The rich 
dynamical behaviors obtained previously in the literature can be regarded 
as the special cases of our present result or Markovian approximation. Par- 
ticularly, we find that the entanglement can stably distribute among all the 
bipartite subsystems if the qubit and its environment can form a bound state 
and the non-Markovian effect is important. Irrespective of how distributes 
the entanglement, it is found that the identity about the entanglement in 
the whole system can be satisfied at any time, which reveals the profound 
physics of the entanglement dynamics. 

4.2 The model of two qubits in two uncorre- 
lated band-gap reservoirs 

The model we studied here is the same as the one used in chapter 3, i.e. two 
qubits in two separate vacuum reservoirs. Because of the independence of the 
two subsystems, we first consider the single "qubit -|- reservoir" subsystem 
which is governed by the following Hamiltonian [6H] 

H = uoa+a- + ^Wfe4afc + "^{gkcr+ak + h.c), (4.1) 

k k 

in which the notations are the same as Sec. 2.2 

As we have discussed in chapter 2, the initial state of the local system 
l</'(0)) = l~) ® l{0}fc)) where |{0}fc) denotes the vacuum state of reservoir, 
does not evolve with time. While for an initial state |0(O)) = |+) ® |{0}fc), 
governed by the Hamiltonian fl4.ip . its time evolution is given by 

|0(t)) = hit) |+, {0},) + Y. h{t) I-, {!},) , (4.2) 
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where |{l}fc) denotes the reservoir state with only one photon in the k-th 
mode. From the Schrodinger equation, we can get the time evolution of the 
excited probability amplitudes in Eq. (14. 2p 

b{t) + iuJob{t) + / 6(r)/(t-r)rfr = 0, (4.3) 



where the kernel function is /(x) = f^ dujJ{uj)e ^"^^ with JiuS) = J2k Id^l ^i 



u- 



ujk) being the spectral density. 

If we define the normalized collective state with one excitation in the 
reservoir as ll) = [^^ hk{t) |{l}fc)]/K^) ^^^ with no excitation in the reser- 
voir as |6)^ = |{0}fc), then Eq. fj42l) can be recast into [¥)1 178] 

|0(t)) = 6(t)|+)|6>^ + 6(t)|-)|i>^, (4.4) 



where 6(t) = -^1 - |6(t)|^ 

According to the above results, the time evolution of a system con- 
sisted of two such subsystems with the initial state |$(0)) = (a;|— ,— ) + 
/3|+,+))|0)^J6)^^isgivenby 

|$(t)) = a I-, -) |6>^^ |6>^^ + /3 m)), m))^ , (4.5) 

where a and /3 are the coefficients to determine the initial entanglement 
in the system. From p = |$(t))($(t)|, one can obtain the time-dependent 
reduced density matrix of the bipartite subsystem qubitl-qubit2 (gig2) by 
tracing over the reservoir variables. It reads 



Pqiq2V') 



1 p? m ' 








Pa*h{tf \ 





P 
















V 







V P*ah*{tf 










X ) 



(4.6) 



where p = \(3b{t)\ h{tY and x = 1 — \(3\'^\b{t)\^ — 2p. By the similar proce- 
dure, it is not difficult to obtain the corresponding reduced density matri- 
ces for other bipartite subsystems like reservoir l-reservoir2 (rir2) and qubit- 
reservoir {qiVi, qir2, q2ri, 52^2)- 
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4.3 Entanglement distribution among bipar- 
tite subsystems 

Follow chapter 3, we use the concurrence [67] to quantify entanglement. 
The concurrence for each bipartite partition can be calculated as Cm = 
niax{0, Qm}, where m denotes the different bipartite partitions and Qm read 



Q,,,, = 2|a/3||6(t)p-2|/36(t)|%)^ 


(4.7) 


Qr,r,=2\amtY-2mt)%itf, 


(4.8) 


Q,^,^ = 2 P\'b{t)m{z = 1,2), 


(4.9) 


Q,,r, = 2\aPbit)\b{t) - 2 m)fKty = Qg,r^- 


(4.10) 


ard to verify that the quantities Qm in Eqs. i 


4.7)-f4.10) 



satisfy an identity 



^qi<12 ~r ^riT2 "T ^ 



Q,,r,-2Q,,r, = 2\aP\, (4.11) 



where 2|a/3| is just the initial entanglement. Eq. (14. lip has been obtained 
in a double J-C model [39], in which the decoherence is absent since each of 
the reservoirs only contains one mode, i.e. J{uj) = g^Siuj — Uq). Surprisingly, 
this identity is still true even in the presence of decoherence. Furthermore, 
one notes that the identity is not dependent of any detail about h{t), which 
is determined by Eq. (14. 3p . This result shows clearly the invariant nature 
of the entanglement. In the following we explicitly discuss the distribution 
behavior of the entanglement by taking the reservoir as a photonic band gap 
(PBG) medium [251 137] and compare it with the previous results. 

For the PBG medium, the dispersion relation near the upper band-edge 
is given by |26j 

Uk = ujc + A{k - kof, (4.12) 

where A ^ ^c/k^, Uc is the upper band-edge frequency and fco is the corre- 
sponding characteristic wave vector. In this case, the kernel function has the 
form [73] 

f{t-r) = 7] f ^e-'^'^^'-^Uk, (4.13) 
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Qi'-i0.5 



Figure 4.1: Entanglement distributions and their time evolutions for the case 
of uq < Uc- The parameters used are Uq = 0.1 Wc and t] = 0.2. 



where rj 



Q^§^ ^3 is a dimensionless constant. In solving Eq. (14. 3 p for b{t), 
Eq. f l4.13p is evaluated numerically. Here we do not make an assumption 
that k can be replaced by ko outside of the exponential [36], as also done 
in Refs. [28| [TH [75]. Thus our result is numerically exact. In the following 
discussion we take Uc as the unit of frequency. 

In Figs. 14.11 and 14. 2[ we show the entanglement distributions and their 
time evolutions for two typical cases of uq < oJc and ujq > Uc, which cor- 
respond to the atomic frequency being located at the band gap and at the 
upper band of the PEG medium, respectively [73|. In the both cases the 
initial entanglement in qiq2 begins to transfer to other bipartite partitions 
with time but their explicit evolutions, in particular the long time behaviors, 
are quite different. In the former case, the entanglement could distribute 
stably among all bipartite partitions. Fig. 14.1( a) shows that after some 
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Qi'-i0.5 



Figure 4.2: Entanglement distributions and their time evolutions for the case 
of uq > Uc- The parameters used are Uq = lO.Ouc and r] = 0.2. 

oscillations, a sizeable entanglement of qiq2 is preserved for the parameter 
regime of 0.3 < a < 1 in the long-time limit. Remarkably, the entanglement 
in qiri{i = 1,2) forms quickly in the full range of a [Fig. 14.1( c)] and domi- 
nates the distribution. On the contrary, only slight entanglement of rir2 is 
formed in a very narrow parameter regime 0.6 < a < 1, as shown in Fig. 
14.1( b). However, when coq is located at the upper band of the PEG medium, 
the initial entanglement in qiq2 is transferred completely to the rir2 in the 
long-time limit, as shown in Fig. 14.21 At the initial stage, qiri{i = 1,2) 
and qir2{q2^i) are entangled transiently, but there is no stable entanglement 
distribution. This result is consistent with that in Refs. [lOl [761 EZ]- It 
is not difficult to understand these results according to Eqs. fl4.7p -( H?T0|) . 
From these equations, one can clearly see that the detailed behavior of the 
entanglement dynamics and its distributions in the bipartite partitions are 
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completely determined by the time-dependent factor |6(t)P of single-qubit 
excited-state population. Fig. 14.31 shows its time evolutions for the corre- 
sponding parameter regimes presented above. We notice that |6(oo)p 7^ 
when cjo is located at the band gap, which means that there is some excited- 
state population in the long-time limit. This is just the population trapping 
which we have discussed in above chapters. Such population trapping just 
manifests the formation of bound states between qi and r^ [5S] . Consequently, 
Qi and Tj are so correlated in the bound states that the initial entanglement 
in qiq2 cannot be fully transferred to rir2. The oscillation during the evo- 
lution is just the manifestation of the strong non-Markovian effect induced 
by the reservoirs. On the contrary, if ojq is located in the upper band, then 
|6(oo)p = and the qubits decay completely to their ground states. In this 
case the bound states between q^ and Tj are absent and the initial entangle- 
ment in qiq2 is completely transferred to the rir2, as clearly shown in Eq. 

dH). 
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Figure 4.3: Time evolution of time-dependent factor of the excited-state 
population for two parameter regimes uq = O.luc (solid line) and lO.Ouc 
(dashed line), rj is taken as 0.2. 
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In addition, in Refs. [10l[76l[77] it was emphasized that the ESD of qiq2 is 
always accompanied with the ESB of rir2- However, this is not always true. 
To clarify this, we examine the condition to obtain the ESD of the qubits 
and the companying ESB of the reservoirs. From Eqs. (14. 7p and (14. 8 p it is 
obvious that the condition is Qqiq^it) < and Qrir2{'t') > at any t and f, 
which means [73] 

\b{t')f < \a\/^/rr\^ < 1 - mf . (4.14) 

In the case without bound states, |6(c)o)| = 0. The condition (I4.14p can 
be satisfied when a < 1/v^. So one can always expect the ESD of the 
qubits and the companying ESB of the reservoirs in the region \a\ < 1/v^, 
as shown in Fig. 14.21 and Refs. [iQl [761 [ZZ] • However, when the bound states 
are available, the situation changes. In particular, when \b(t)\ > | in the 
full range of time evolution, no region of a can make the condition (I4.14p to 
be satisfied anymore. For clarification, we present three typical behaviors of 
the entanglement distribution in Fig. 14.41 In all these cases the bound states 
are available. Fig. 14.41 (a) shows the situation where the entanglement is 
stably distributed among all of the bipartite subsystems. In Fig. 14.41 (b) the 
entanglement of rir2 shows ESB and revival. However, the entanglement in 
qiq2 does not exhibit ESD. This is the example that the ESD in qiq2 is not 
accompanied with the ESB in rir2. Fig. 14.4( c) shows another example that 
while the entanglement of qiq2 shows ESD and revival [29], the entanglement 
of rir2 does not show ESB but remains to be zero. 

The above discussion is general and is not dependent of the explicit form 
of the reservoir. To confirm this, we consider the radiation field in free space. 
The spectral density has the Ohmic form J{uj) = r]u exp{—uj / A) , which can 
be obtained from the free-space dispersion relation u = ck. One can verify 
that the condition for the formation of bound states is: Uq — r/A < [66]. 
In Fig. 14.51 we plot the results in this situation. The previous results can 
be recovered when the bound states are absent. On the contrary, when the 
bound states are available, a stable entanglement is established among all 
the bipartite partitions. Therefore, we argue that the stable entanglement 
distribution resulted from the bound states is a general phenomenon in open 
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Figure 4.4: (Color online) Entanglement evolution when a = l/v2 (a), 
a = 0.57 (b), and a = 0.28 (c). The parameters used here are the same as 
Fie. SH 
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quantum system when the non-Markovian effect is talcen into account. 

4.4 Summary 

In summary, we have studied the entanglement distribution among all the 
bipartite subsystems of two qubits embedded into two independent reser- 
voirs. It is found that the entanglement can be stably distributed in all the 
bipartite subsystems and they satisfy an identity about the entanglement. 
This identity is shown to be independent of any detail of the reservoir and 
its coupling with the qubit, which affect only the explicit time evolution be- 
havior and the final distribution. The result shows the physical nature of the 
entanglement and has a significant implication for the quantum information 
processing. 
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Chapter 5 

Summary and outlook 



In this thesis we studied the decoherence dynamics of open quantum system. 
We found that the decoherence would be greatly suppressed if the bound 
state is formed under the non-Markovian dynamics. 

We model our system as two-level atoms in vacuum reservoirs. After 
numerically solving the coupled equations, we studied the non-Markovian 
effect on the decoherence dynamics of the quantum system. Compared with 
the results obtained under the Markovian approximation, we found that the 
environment has two effects on the quantum system of interest: dissipation 
effect, which degrades the quantum coherence, and backaction effect, which 
compensates the quantum system of lost coherence. The competition of these 
two effects results in the rich decoherence dynamic behaviors. Our results 
show explicitly, in the weak coupling regime, the widely used Markovian 
approximation is applicable, while, in the strong coupling regime, the non- 
Markovian effect makes the dynamical process oscillate for a certain time 
which is a manifestation of information or energy flowing back and forth be- 
tween quantum system and the environment. We also find that the quantum 
coherence can be preserved in the long time limit. Physically, it is attribute 
to the formation of a bound state between the quantum system and its local 
reservoir. We give explicitly a condition to judge the formation of bound 
state for any kind of vacuum reservoir including the widely studied PBG 
medium. 

Due to the preservation of quantum coherence of the single system, we 
have revealed that the entanglement for a composite two-qubit system can 
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also be preserved in the steady state. We give explicitly the mechanism of 
entanglement preservation, i.e. the fulfillment of non-Markovian effect and 
formation of bound states. The mechanism we given can explain the results 
obtained in the previous works. When the non-Markovian effect is neglected, 
the phenomenon of ESD of the qubits is reproduced. The phenomenon of 
ESD and its revival can be obtained when the non-Markovian effect is con- 
tained while the bound state is not available. In particular, the entanglement 
preservation when the atoms are placed in the PBG mediums reported in 
Ref. [28] can be explained as the fulfillment of the above two conditions. In 
a word, we have given a clear clue on how to preserve entanglement in the 
steady state. 

Considering the environment as a whole, we also investigated the entan- 
glement distribution among all the bipartite subsystems. We found that the 
entanglement can be stably distributed among all bipartite partitions of the 
whole system when the bound states are formed. It is particularly interesting 
to find that the entanglement in different bipartite partitions always satisfies 
an identity. This identity is independent on the explicit dynamics process. 
Our unified treatment includes the previous results in the literature as special 
cases. When the bound state is absent and the Markovian approximation is 
applicable, the result reported in Ref. [IQ] that the entanglement transfer 
from the qubits to the reservoirs is recovered. Our work give a thorough un- 
derstanding of entanglement distribution among quantum systems and their 
environments. 

There are many open issues relevant to the subject of this thesis. For 
example, the mechanism of formation of bound state is still unclear, we think 
it is a kind of quantum phase transition when the bound state is formed. How 
to reveal the relationship between the bound state and the quantum phase 
transition is an open question. What's more, in this thesis we have shown 
that the non-Markovian effect can rescue the entanglement. Does the non- 
Markovian effect also can rescue certain missions of quantum information 
processing, for example, quantum teleportation and quantum dense coding, 
in noisy quantum channels? This is still an open question. 



46 



Bibliography 



[1] H.-P. Breuer and F. Petruccione, The Theory of Open Quantum Systems 
(Oxford University Press, Oxford, 2002). 

[2] M. A. Nielsen and I. L. Chuang, Quantum, Computation and Quantum Infor- 
mation (Cambridge University Press, Cambridge, 2000). 

[3] A. Einstein, B. Podolsky, and N. Rosen, Can Quantum- Mechanical Descrip- 
tion of Physical Reality Be Considered Complete?, Phys. Rev. 47, 777(1935). 

[4] J. S. Bell, On the Einstein-Podolsky-Rosen Paradox, Physics (NY) 1, 195 
(1964). 

[5] W. H. Zurek, Decoherence, einselection, and the quantum origins of the clas- 
sical , Rev. Mod. Phys. 75, 715 (2003). 

[6] C. H. Bennett, F. Bessette, G. Brassard, L. Salvail, and J. Smolin, Experi- 
mental quantum cryptograph, J. Cryptol. 5, 3 (1992). 

[7] C. H. Bennett, G. Brassard, C. Crepeau, R. Jozsa, A. Peres, and W. K. Woot- 
ters, Teleporting an unknown quantum state via dual classical and Einstein- 
Podolsky-Rosen channels, Phys. Rev. Lett. 70, 1895 (1993). 

[8] J.-W. Pan, D. Bouwmeester, H. Weinfurter, and A. Zeilinger, Experimental 
Entanglement Swapping: Entangling Photons That Never Interacted, Phys. 
Rev. Lett. 80, 3891 (1998). 

[9] P. W. Shor, in Proceedings of the 35th Annual Symposium on the Founda- 
tions of Computer Science, Santa Fe, 1994, edited by S. Goldwasser (IEEE 
Computer Society Press, Los Alamo, CA, 1994). 
[10] L. K. Grover, Quantum Mechanics Helps in Searching for a Needle in a 

Haystack, Phys. Rev. Lett. 79, 325 (1997). 
[11] A. K. Ekert, Quantum cryptography based on Bell's theorem, Phys. Rev. Lett. 

67, 661 (1993). 
[12] C. W. Gardiner and P. ZoUer, Quantum Noise (Springer, Berlin, 2004). 
[13] K. Zyczkowski, P. Horodecki, M. Horodecki, and R. Horodecki, Dynamics of 

quantum entanglement, Phys. Rev. A 65, 012101 (20001). 
[14] T. Yu and J. H. Eberly, Finite-Time Disentanglement Via Spontaneous Emis- 
sion, Phys. Rev. Lett. 93, 140404 (2004); Science 323, 598 (2009). 
[15] M. Ikram, Fu-li Li, and M. S. Zubairy, Disentanglement in a two-qubit system 
subjected to dissipation environments, Phys. Rev. A 75, 062336 (2007). 



47 



[16] M. P. Almeida, F. de Melo, M. Hor-Meyll, A. Salles, S. P. Walborn, P. H. S. 
Ribeiro, and L. Davidovich, Experimental observation of environment-induced 
sudden death of entanglement, Science 316, 579 (2007). 

[17] J. Laurat, K. S. Choi, H. Deng, C. W. Chou, H. J. Kimble, Heralded Entan- 
glement between Atomic Ensembles: Preparation, Decoherence, and Scaling, 
Phys. Rev. Lett. 99, 180504 (2007). 

[18] M. Fujita, S. Takahashi, Y. Tanaka, T. Asano, and S. Noda, Simultaneous In- 
hibition and Redistribution of Spontaneous Light Emission in Photonic Crys- 
tals, Science, 308, 1296 (2005). 

[19] F. Dubin, D. Rotter, M. Mukherjee, C. Russo, J. Eschner, and R. Blatt, 
Photon Correlation versus Interference of Single- Atom Fluorescence in a Half- 
Cavity, Phys. Rev. Lett. 98, 183003 (2007). 

[20] F. H. L. Koppens, D. Klauser, W. A. Coish, K. C. Nowack, L. P. Kouwen- 
hoven, D. Loss, and L. M. K. Vandersypen, Universal Phase Shift and Non- 
exponential Decay of Driven Single-Spin Oscillations, Phys. Rev. Lett. 99, 
106803 (2007). 

[21] D. Mogilevtsev, A. P. Nisovtsev, S. Kilin, S. B. Cavalcanti, H. S. Brandi, and 
L. E. Oliveira, Driving-Dependent Damping of Rabi Oscillations in Two-Level 
Semiconductor Systems, Phys. Rev. Lett. 100, 017401 (2008). 

[22] J.-S. Xu, C.-F. Li, M. Gong, X.-B. Zou, C.-H. Shi, G. Chen, and G.-C. Guo, 
Experimental Demonstration of Photonic Entanglement Collapse and Revival, 
Phys. Rev. Lett. 104, 100502 (2010). 

[23] M. Scala, B. Militello, A. Messina, S. Maniscalco, J. Piilo, and K.-A. Suomi- 
nen, Population trapping due to cavity losses, Phys. Rev. A 77, 043827 (2008). 

[24] P. Kaer, T. R. Nielsen, P. Lodahl, A. -P. Jauho, and J. M0rk, Non-Markovian 
Model of Photon- Assisted Dephasing by Electron- Phonon Interactions in a 
Coupled Quantum-DotCCavity System, Phys. Rev. Lett. 104, 157401 (2010). 

[25] E. Yablonovitch, Inhibited Spontaneous Emission in Solid-State Physics and 
Electronics, Phys. Rev. Lett. 58, 2059 (1987). 

[26] S. John and T. Quang, Spontaneous emission near the edge of a photonic 
band gap, Phys. Rev. A 50, 1764 (1994). 

[27] P. Lambropoulos, G. M. Nikolopoulos, T. R. Nielsen, and S. Bay, Fundamental 
quantum optics in structured reservoirs. Rep. Prog. Phys. 63, 455 (2000). 

[28] B. Bellomo, R. LoFranco, S. Maniscalco, and G. Compagno, Entanglement 
trapping in structured environments, Phys. Rev. A 78, 060302(R) (2008). 

[29] B. Bellomo, R. LoFranco, and G. Compagno, Non-Markovian Effects on the 
Dynamics of Entanglement, Phys. Rev. Lett. 99, 160502 (2007); Entangle- 
ment dynamics of two independent qubits in environments with and without 
memory, Phys. Rev. A 77, 032342 (2008). 

[30] K.-L. Liu and H.-S. Goan, Non-Markovian entanglement dynamics of quan- 
tum continuous variable systems in thermal environments, Phys. Rev. A 76, 
022312 (2007). 



48 



[31] J.-H. An and W.-M. Zhang, Non-Markovian entanglement dynamics of noisy 

continuous variable quantum channels , Phys. Rev. A 76, 042127 (2007). 
[32] J.-H. An, M. Feng, and W.-M. Zhang, Non-Markovian decoherence dynamics 

of entangled coherent states , Quant. Inf. Comp. 9, 0317 (2009). 
[33] J.-H. An, Y. Yeo, and C. H. Oh, Exact decoherence dynamics of a single-mode 

optical field, Ann. Phys. (N.Y.) 324, 1737 (2009). 
[34] J.-H. An, Y. Yeo, W.-M. Zhang, and C. H. Oh, Entanglement oscillation 

and survival induced by non-Markovian decoherence dynamics of entangled 

squeezed-state, J. Phys. A: Math. Theor. 42, 015302 (2009). 
[35] S. John and J. Wang, Quantum electrodynamics near a photonic band gap: 

Photon bound states and dressed atoms, Phys. Rev. Lett. 64, 2418 (1990). 
[36] M. Woldeyohannes and S. John, Coherent control of spontaneous emission 

near a photonic band edge: A qubit for quantum computation, Phys. Rev. A 

60, 5046 (1999). 
[37] P. Lodahl, A. F. van Driel, I. S. Nikolaev, A. Irman, K. Overgaag, D. Van- 

niaekelbergh, and W. L. Vos, Controlling the dynamics of spontaneous emis- 
sion from quantum dots by photonic crystals. Nature (London) 430 , 654 

(2004). 
[38] M. Miyamoto, Bound-state eigenenergy outside and inside the continuum for 

unstable multilevel systems, Phys. Rev. A 72, 063405 (2005). 
[39] M. Yonag, T. Yu, and J. H. Eberly, Pairwise Concurrence Dynamics: A 

Four-Qubit Model, J. Phys. B. 40, S45 (2007). 
[40] C. E. Lopez, G. Romero, F. Lastra, E. Solano, and J. C. Retamal, Sudden 

Birth versus Sudden Death of Entanglement in Multipartite Systems, Phys. 

Rev. Lett. 101, 080503 (2008). 
[41] S. Mancini and R. Bonifacio, Quantum Zeno-like effect due to competing de- 
coherence mechanisms, Phys. Rev. A 64, 042111 (2001). 
[42] L. M. Duan and G. C. Guo, Preserving Coherence in Quantum Computation 

by Pairing Quantum Bits, Phy. Rev. Lett. 79, 1953 (1997). 
[43] L. Viola and S. Lloyd, Dynamical suppression of decoherence in two-state 

quantum systems, Phys. Rev. A 58, 2733 (1998). 
[44] D. Vitali and P. Tombesi, Using parity kicks for decoherence control, Phys. 

Rev. A 59, 4178 (1999). 
[45] P. Facchi, S. Tasaki, S. Pascazio, H. Nakazato, A. Tokuse, and D. A. Lidar, 

Control of decoherence: Analysis and comparison of three different strategies, 

Phys. Rev. A 71, 022302 (2005). 
[46] B. M. Garraway, Nonperturbative decay of an atomic system in a cavity, Phys. 

Rev. A 55, 2290 (1997). 
[47] S. Maniscalco and F. Petruccione, Non-Markovian dynamics of a qubit, Phys. 

Rev. A 73, 012111 (2006). 
[48] J. Piilo, S. Maniscalco, K. Harkonen, and K.-A. Suominen, Non-Markovian 

Quantum Jumps, Phys. Rev. Lett. 100, 180402 (2008). 



49 



[49] W. Koch, F. Grofimann, J. T. Stockburger, and J. Ankerhold, Non-Markovian 

Dissipative Semiclassical Dynamics, Phys. Rev. Lett. 100, 230402 (2008). 
[50] H.-P. Breuer and B. Vacchini, Quantum Semi-Markov Processes, Phys. Rev. 

Lett. 101, 140402 (2008). 
[51] P. Rebentrost, L Serban, T. Schulte-Herbriiggen, and F. K. Wilhclm, Optimal 

Control of a Qubit Coupled to a Non-Markovian Environment, Phys. Rev. 

Lett. 102, 090401 (2009). 
[52] H.-P. Breuer, E.-M. Laine, and J. Piilo, Measure for the Degree of Non- 
Markovian Behavior of Quantum Processes in Open Systems, Phys. Rev. Lett. 

103, 210401 (2009). 
[53] D. Chruscinski and A. Kossakowski, Non-Markovian Quantum Dynamics: 

Local versus Nonlocal, Phys. Rev. Lett. 104, 070406 (2010). 
[54] D. P. DiVincenzo and D. Loss, Rigorous Born approximation and beyond for 

the spin-boson model, Phys. Rev. B 71, 035318 (2005). 
[55] C. J. Myatt, B. E. King, Q. A. Turchette, C. A. Sackett, D. Kielpinski, W. M. 

Itano, C. Monroe, and D. J. Wineland, Decoherence of quantum superpositions 

through coupling to engineered reservoirs. Nature (London) 403, 269 (2000). 
[56] Q. A. Turchette, C. J. Myatt, B. E. King, C. A. Sackett, D. Kielpinski, W. 

M. Itano, C. Monroe, and D. J. Wineland, it Decoherence and decay of 

motional quantum states of a trapped atom coupled to engineered reservoirs, 

Phys. Rev. A 62, 053807 (2000). 
[57] S. Diehl, A. Micheh, A. Kantian, B. Kraus, H. P. Bchler, and P. Zoller, 

Quantum states and phases in driven open quantum systems with cold atoms, 

Nat. Phys. 4, 878 (2008). 
[58] A. J. Leggett, S. Chakravarty, A. T. Dorsey, M. P. A. Fisher, A. Gary, and 

W. Zwerger, Dynamics of the dissipative two-state system. Rev. Mod. Phys. 

59, 1 (1987). 
[59] Q.-J. Tong, J.-H. An, H.-G. Luo, and C. H. Oh, Decoherence suppression of 

a dissipative qubit by non-Markovian effect, J. Phys. B: At. Mol. Opt. Phys. 

43 155501 (2010). 
[60] B. L. Hu, J. P. Paz, and Y. Zhang, Quantum Brownian motion in a general 

environment: Exact master equation with nonlocal dissipation and colored 

noise, Phys. Rev. D 45, 2843 (1992). 
[61] U. Weiss, Quantum Dissipative Systems, 3nd ed. (World Scientific, Singapore, 

2008). 
[62] J. P. Paz and A. J. Roncaglia, Dynamical phases for the evolution of the 

entanglement between two oscillators coupled to the same environment, Phys. 

Rev. A 79, 032102 (2009). 
[63] I. E. Linington and B. M. Garraway, Dissipation control in cavity QED with 

oscillating mode structures, Phys. Rev. A 77, 033831 (2008). 
[64] D. Porras, F. Marquardt, J. von Delft, and J. L Cirac, Mesoscopic spin-boson 

models of trapped ions, Phys. Rev. A 78, OIOIOI(R) (2008). 



50 



[65] I. Buluta and F. Nori, Quantum Simulators, Science 326, 108 (2009). 

[66] Q.-J. Tong, J.-H. An, H.-G. Luo, and C. H. Oh, Mechanism, of entanglem,ent 
preservation, Phys. Rev. A 81, 052330 (2010). 

[67] W. K. Wootters, Entanglement of Formation of an Arbitrary State of Two 
Qubits, Phys. Rev. Lett. 80, 2245 (1998). 

[68] M. O. Scuhy and M. S. Zubairy, Quantum Optics (Cambridge University 
Press, Cambridge, 1997). 

[69] A. Walhaff, D. I. Schuster, A. Blais, L. Frunzio, R.- S. Huang, J. Majer, S. 
Kumar, S. M. Girvin, and R. J. Schoelkopf, Strong coupling of a single pho- 
ton to a superconducting qubit using circuit quantum electrodynamics, Nature 
(London) 431, 162 (2004). 

[70] K. Hennessy, A. Badolato, M. Winger, D. Gerace, M. Atatiire, S. Guide, S. 
Fait, E. L. Hu, and A. Imamoglu, Quantum nature of a strongly coupled single 
quantum dotCcavity system. Nature (London) 445, 896 (2007); A. Faraon, A. 
Majumdar, H. Kim, P. Petroff, and J. Vuckovic, Fast Electrical Control of 
a Quantum Dot Strongly Coupled to a Photonic- Crystal Cavity, Phys. Rev. 
Lett. 104, 047402 (2010). 

[71] Q. A. Turchette, C. J. Myatt, B. E. King, C. A. Sackett, D. Kielpinski, W. M. 
Itano, C. Monroe, and D. J. Wineland, Decoherence and decay of motional 
quantum states of a trapped atom coupled to engineered reservoirs, Phys. Rev. 
A 62, 053807 (2000); S. Diehl, A. Micheh, A. Kantian, B. Kraus, H. P. 
Biichler, and P. Zoller, Quantum states and phases in driven open quantum 
systems with cold atoms, Nat. Phys. 4, 878 (2008). 

[72] I. E. Linington and B. M. Garraway, Control of atomic decay rates via ma- 
nipulation of reservoir mode frequencies, J. Phys. B: At. Mol. Opt. Phys. 39, 
3383 (2006). 

[73] Q.-J. Tong, J.-H. An, H.-G. Luo, and C. H. Oh, Entanglement distribution 
and its invariance, I arXiv: 1005.1001] submitted. 

[74] M. Al-Amri, G.-x. Li, R. Tan, and M. S. Zubairy, Sudden death and birth of 
entanglement in photonic crystals, Phys. Rev. A 80, 022314 (2009). 

[75] F.-Q. Wang, Z.-M. Zhang, and R.-S. Liang, Finite-time decoherence could be 
suppressed efficiently in a photonic crystal, Phys. Rev. A 78, 042320 (2008). 

[76] J. Zhou, C.-J. Wu, M.-Y. Zhu, and H. Guo, Non-Markovian dynamics of 
entanglement for multipartite systems, J. Phys. B: At. Mol. Opt. Phys. 42, 
215505 (2009). 

[77] Z. Y. Xu and M. Feng, Sudden death and birth of entanglement beyond the 
Markovian approximation, Phys. Lett. A 373, 1906 (2009). 



51 



Publication list 



1. Qing-Jun Tong, Jun-Hong An, Hong-Gong Luo, and C. H. Oh, Mechanism 
of entanglement preservation, Phys. Rev. A 81, 052330 (2010). 

2. Juan-Juan Chen, Jun-Hong An, Qing-Jun Tong, Hong-Gang Luo, and C. 
H. Oh, Non-Markovian effect on the geometric phase of a dissipative qubit, 
Phys. Rev. A 81, 022120 (2010). 

3. Qing-Jun Tong, Jun-Hong An, Hong-Gong Luo, and C. H. Oh, Decoher- 
ence suppression of a dissipative qubit by non-Markovian effect, J. Phys. B: 
At. Mol. Opt. Phys. 43 155501 (2010). 

4. Qing-Jun Tong, Jun-Hong An, Hong-Gang Luo, and C. H. Oh, Entangle- 
ment distribution and its invariance, .arXiv:1005.1001^ submitted. 



52 



Acknowledgments 



As this thesis draws to a close, I give my great thanks to my advisor, 
Professor Jun-Hong An for his patient instruction and warm-hearted helps 
during the past years. Without his guidance and help, it is impossible for 
me to go into the study of decoherence dynamics of open quantum system 
and accomplish this thesis. 

Furthermore, I express my great thanks to Professor Hong-Gang Luo for 
his illuminating discussions and instructive helps. I also acknowledge Dr. 
Zheng-Xiao Zhou, Dr. Lin Li, Juan- Juan Chen and many other colleagues 
for their valuable helps and discussions during my research works. 

Finally, I can not forget the long-term support of my parents, my brother 
and my lover who have shown patience and understanding to the entire period 
of my work. 



53 



